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ABSTRACT 

Context. Understanding the source of systematic errors in photometry is essential for their calibration. 

Aims. We investigate how photometry performed on difference images can be influenced by errors in the photometric scale factor. 
Methods. We explore the equations for difference image analysis (DIA) and we derive an expression describing how errors in the 
difference flux, the photometric scale factor and the reference flux are propagated to the object photometry. 

Results. We find that the error in the photometric scale factor is important, and while a few studies have shown that it can be at a 
significant level, it is currently neglected by the vast majority of photometric surveys employing DIA. 

Conclusions. Minimising the error in the photometric scale factor, or compensating for it in a post-calibration model, is crucial for 
reducing the systematic errors in DIA photometry. 
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1. Introduction 


2. Equations 


The technique of difference image ana l 


Alard & Lur)tonlll998 ; Alardll2000t lBramichll2008l: IBecker et al.l 


(DIA; 


2Q12t iBramich et al.l 1201 3h is based on matching a reference 

image to a target image by modelling the differences in align¬ 
ment, point-spread function (PSF), exposure time, atmospheric 
extinction, and sky background between them. Specifically, a 
convolution kernel is used to model the first four differences (to 
within a small translational offset) while an additive differential 
background is used to model the last. The reference image is 
transformed by convolving it with the kernel and adding the dif¬ 
ferential background, and the result is subtracted from the target 
image to create a difference image. All non-varying sources are 
fully subtracted in the difference image, leaving signal only for 
sources that have varied in brightness (or possibly position). 

DIA is being increasingly used for precision photometry 
and transient detecti on in a wide range of photo metric surveys 
(e.g. P anSTARRS - iKaiser et al.ll2nn^ OGEE - lUdalskil 120031 


LSST -llvezic eyal.l2008l RoboNet-II - iTsanras et al.l20ld OES 
- lAlsubai et al.l ^131 etc.). Furthermore, these surveys are in¬ 
vesting substantial efforts into post-calibration in order to min¬ 
imise the systematic noise in the survey photometry which af¬ 
fects important aspects suc h as the detection limits, homogene¬ 
ity and comple teness (e.g. IStubbs et al.ll201()[ lOfek et aTll2012i 
IWittman et al.ll201^ etc.). Therefore, it is crucial to understand 
how the difference images, on which the photometry is per¬ 
formed, are created and how systematic errors in the difference 
images themselves can affect the photometry. However, there is 
no study in the literature on the systematic errors specific to DIA. 
In this research note, we have opened the investigation into sys¬ 
tematic errors in DIA by exploring the effect that an error in the 
kernel sum, known as the photometric scale factor, may have on 
the photometry. 


The target image I is modelled as the convolution of a refer¬ 
ence image R with a convolution kernel plus a differential back¬ 
ground. If one considers the kernel as the product of a photomet¬ 
ric scale factor P and a normalised (by its sum) kernel K, then 
the model target-image M is defined as follows: 

M ^ P{R®K) + B (1) 

The corresponding difference image D is the image of model 
residuals given by: 

D^I-M 

-P{R®K)-B 

Now assume that an object, consisting of a source of interest 
and a blend, has a true flux /(f) on the photometric scale of the 
reference image that is given by: 

where k{t) is a function of time t that represents any variability 
in the source and kb is the true blend ratio. Without loss of gen¬ 
erality, adopt k{0) = 0 so that the quantity fs represents the true 
source flux at f = 0. The true blend flux kb fs makes the source 
appear brighter by a factor of (1 -i- kb) and reduces the apparent 
fractional flux-amplitude of any source variability by the same 
factor. 

It follows that the true object flux /* on the reference image 
taken at f = 0 is; 


/«-(!+ h) fs 


(4) 
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and that the true object flux //(f) on the target image taken at 
time t is: 

//(f) = Pf(t) (5) 

By considering equation|2] the true object flux foil) on the dif¬ 
ference image may be computed as: 

fait) ^ fi(t) - P fR (6) 

By substituting equation |5] into equation |6] and rearranging, /(f) 
may be written in terms of the measurable quantities /r and foil) 
as: 


fit)^fR+fD(t)IP 


(7) 


This is the equation used to con vert difference fluxes to total 
fluxes and thence magnitudes (e.g. lBramich et al1l201 Ih . 

Unfortunately, the modelling of the target image is never per¬ 
fect and the difference image consequently suffers from small 
systematic errors that are propagated to the photometry. An er¬ 
ror in the fitted differential background which results in a non¬ 
zero background in the difference image is trivially accounted 
for, at the expense of a little extra variance in the photometry, 
by including the background as a parameter in the method used 
for performing the photometry on the difference image (i.e. PSF 
fitting or aperture photometry). However, an error in the PSF 
matching (both shape and scale) produces systematic residuals 
at the object positions that are more difficult to mitigate at the 
image processing stage. Aperture photometry is agnostic to mis¬ 
matches in PSF shape, but will be affected by an error in the fit¬ 
ted photometric scale factor. PSF photometry on the other hand 
is sensitive to mismatches in both PSF shape and scale but will 
provide photometry with smaller variance than aperture photom¬ 
etry when stochastic noise dominates. For the above reasons, we 
ignore the error in the fitted differential background, we treat 
the error in the fitted normalised kernel as part of the error in¬ 
troduced into the photometry by the measurement process on 
the difference image, and we assume that the fitted photometric 
scale factor suffers from a small fractional error. 

We use P', K', and B' to represent the fitted photometric 
scale factor, kernel, and differential background, respectively. 
Employed in equation |2] these give us the difference image D' 
with systematic errors. If we assume that the difference flux is 
measured with an error of sd fi(t) due to stochastic noise and/or 
the error in K', and if we also assume that the error in B' can be 
successfully accounted for, then we obtain the following mea¬ 
sured difference flux on D' for the object: 

fD'it)-Mt)-p'fR+SDm ( 8 ) 

Adopting the expression P' - {\-i-ep)P for the relation between 
P' and P, and using equations[2to|5] we may derive: 

//^'(0 = F (( 1 -H /j/j) k(t) + {sd - fip) ( 1 h- kb)) fs (9) 


If we now assume that the method of performing photometry on 
the reference image yields a fractional flux error of sr in the ref¬ 
erence flux, different to sd because of the different nature of the 
reference image and/or method used, then the measured refer¬ 
ence flux is /r - (1-1- er) fR. Using //, //,(/ and P' in equa¬ 
tion |7] and doing some algebra, yields the following expression 
for the measured object flux /'(f) on the photometric scale of the 
reference image: 


/'(f)= 1 + 


k{f) W1 + So 

{l+d){l+kb)j\l+Ep 


(l+d)il+kb)fs (10) 
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where: 


6 ^ Er 


1 -H £p 

i + Ed 


( 11 ) 


Equation [TO] describes how the measured flux of a constant 
source (k(f) = 0 for all f), or a variable source {k{t) + 0 for 
at least some f), is distorted by the errors ed, Ep and er. The 
equation has been written in the form above to facilitate direct 
comparison to equation |3] representing the true object flux. The 
ratio of the measured to the true object flux is: 


fit) 

m 


1 H- 


l + iw 

\+h 


\ + Ed 
1 -I- Ep 


( 12 ) 


In magnitudes, eauations[3]and[T0lbecome: 


m(f) = -2.51og(/s)-2.51og 1 + 


k(t) \ 
l+kbj 


-2.5log(l+kb) (13) 


».'(r) = - 2.5M/,)-2.5log(1 + 

- 2.5 log(l -H S//) H-2.5 log(l -H sp) 

-2.51og(l -Hd)-2.51og(l H-ki) 

where m(t) and m'{t) are the true and measured object magni¬ 
tudes, respectively, on the magnitude scale of the reference im¬ 
age. 


3. Discussion 

The difference flux ///(f) is a quantity that is measured for each 
object on each difference image, and therefore ed is specific to 
the object and difference image under consideration. However, 
for multiple difference images, any systematic (as opposed to 
stochastic) component in ed that is a function of either an object 
property (e.g. colour) and/or an image property (e.g. pixel coor¬ 
dinates) may be estimated by solving for the appropriate mag¬ 
nitude offsets using the DIA photometry of all of the constant 
objects in the corresponding ta rget images . This i s the approach, 
develo ped by authors such as iHonevcut^ (Il992h and iManfroidI 
(Il995h . that is starting to be adopted by many surveys as the 
standard procedure for performing a post-calib ration of the pho¬ 
tometric data (e.g. iPadmanabhan et akl l2008h . In this respect, 
post-calibration of DIA photometry is no different than the post¬ 
calibration of photometry performed directly on the target im¬ 
ages. The appropriate magnitude offsets to be determined from 
the constant objects are represented by the term -2.5 log(l -H ed) 
in equation [14] and their absolute values are usually of the order 
of ~ 1 - 30 mmag. 

The reference flux /« is a quantity that is measured for each 
object on the reference image. Therefore the error er is inde¬ 
pendent of the target image (or time) and it affects the photome¬ 
try of constant objects by making them systematically too bright 
(er > 0) or too faint (er < 0). Variable objects suffer this same 
systematic error and in addition their fractional flux-amplitude 
of variation is either systematically amplified {er < 0) or re¬ 
duced (er > 0) by a factor of(lH-d)~(lH- er) to first order. In 
this respect, the effect of er is equivalent to that of an extra blend 
flux. The absolute value of er is usually in the range of typical 
photometric precisions of ~ 0.1 - 5%. Unless the DIA pho¬ 
tometry has a very small stochastic noise component (^0.01%), 
Er is indistinguishable from kb since its effect on /'(f) only dif¬ 
fers from that of kb to second order (see equations 
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We do not consider any further the intricacies of disentangling 
the source flux, blend ratio and reference flux error. Simply we 
note that when the blend ratio can be determined (e.g. by using 
external information, or for certain types of variability such as 
microlensing events), the estimated blend ratio is related to 
the true blend ratio by ~ (1 + er) kj, + er. 

When fitting the model target-image, the photometric scale 
factor P is typically assumed to be spatially-invariant and hence 
characterised by a single number, althoug h P may also be mod - 
elled as a function of detector coordinates (iBramich et aLll^OOh . 
Regardless of how P is modelled, the error Ep will be different 
for each target image since it is determined on a per-image basis 
from noisy images. Other effects such as flat-fielding errors, or 
changing non-uniform atmopsheric extinction (i.e. clouds and/or 
airmass gradients), may conspire to make sp specific to the ob¬ 
ject and target image under consideration. As with ed, some of 
the systematic components in Ep may be estimated by solving 
for the appropriate magnitude offsets using the DIA photometry 
of all of the constant objects in the target images. In fact, the 
systematic components in sp that are also common to those in 
Ed may be absorbed into the magnitude offsets represented by 
the term -2.5 log(l -H ed) in equation[T4l However, the remain¬ 
ing systematic components in Ep, if left uncorrected, will cause 
errors in the DIA photometry that have the potential to be misin¬ 
terpreted as real signals. 


4. Summary And Recommendations 


The vast majority of photometric surveys employing DIA do not 
consider R to be a source of error since it is assumed that P has 
been precisely determined, which is equivalent to assuming that 
the differences in atmospheric extinction between the reference 
and target images have been perfectly accounted for by the DIA 
modelling itself. However, in this research note we have shown 
that any errors in P that do exist will have an important impact 
on the DIA photometry. We therefore strongly recommend that 
it becomes standard procedure to assess, and if necessary correct 
for, the effect of the mean error in P for each target image on the 
DIA photometry (encapsulated by the term 2.5 log(l -i- Ep(t)) in 
equation O. 

One method to do this involves htting a post-calibration pho¬ 
tometric model including a set of per-image magnitude offsets 
2.5 log(l -H Ep{t)) to the DIA photometry of all of the constant 
objects in the target images. Note that the post-calibration model 
should also incl ude the magnitudes of the con stant objects as free 
parameters (see IBramich & Freudlindl2012h and any other rele¬ 
vant terms such as -2.5 log(l + eq)- If the variations in the per- 
image magnitude offsets are found to be smaller than the level of 
the stochastic noise in the best object photometry, then they may 
be dropped from the post-calibration model. However, if they are 
deemed to be significant, then the per-image magnitude offsets 
2.5 log(l -H Ep(t)) may be used to correct the DIA photometry of 
all of the objects for the photometric error introduced by Ep(t). 

This technique has started to be used in the series of pa- 
pers on variable stars in globular clusters by t he lead author (e.g. 
iKains et ^120131 lArellano Ferro et ^ 20131 etc, which employ 
the methodology of Bramich & Freudlind (1201 jh i. and in these 
works the relevant per-image magn itude offsets are fo und to be 
at the 0.1-2% level (see Figure 1 of iKains et alJl2015h . This ex¬ 
ample clearly demonstrates that DIA does not always perform 
the photometric matching between images to a precision that is 
below the stochastic noise in the photometric measurements and 
it serves to emphasise how crucial it is to account for the mean 


error in P for each target image in order to minimise the associ¬ 
ated systematic errors. 

Finally we caution that when low-level (^1-2%) suspected 
signals occur in the DIA photometry of an object that cannot be 
confirmed either by the detection of a repeating signal (for peri¬ 
odic signals) or by independent observations taken at the same 
epoch, a careful analysis of the images on which the suspected 
signal was detected is warranted. Checks should be performed 
for the presence of clouds (e.g. light cirrus) that may have caused 
non-uniform atmospheric extinction across the field-of-view that 
also evolves throughout the time-series with the cloud movement 
since this can cause smooth temporal variations in Ep that are 
different for each object and that will manifest themselves as 
smooth variations in the object light curves. Adopting a spatially 
variable photometric scale factor in the model target-image for 
the DIA may partly mitigate this problem. 
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